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ABSTRACT 

Laser anemometry in unseeded flows based on the measurement 
of the spectrum of Rayleigh scattered laser light is reviewed. 
The use of molecular scattering avoids the well-known problems 
(particle lag, biasing effects, seed generation, seed injection) of 
seeded flows. The fundamental limits on velocity measurement 
accuracy are determined using maximum likelihood methods. 
Measurement of the Rayleigh spectrum with scanning 
Fabry— Perot interferometers is analyzed and accuracy limits are 
established for both single pass and multi-pass configurations. 
Multi-pass configurations have much higher selectivity and are 
needed for measurements where there is a large amount of excess 
noise caused by stray laser light. It is shown that Rayleigh 
scattering is particularly useful for supersonic and hypersonic 
flows. The results of the analysis are compared with 
measurements obtained with a Rayleigh scattering diagnostic 
developed for study of the exhaust plume of a small 
hydrogen— oxygen rocket, where the velocities are in the range 
1000 to 5000 m/sec. 


INTRODUCTION 

Conventional laser anemometry is based on the scattering of 
light from particles, usually in the size range 0.5 to 2 
micrometers diameter. In spite of the general success of this 
technique, there are severe limitations in certain applications. 
The introduction of suitable seed material into the flow is often 
difficult, or simply not feasible, such as in the study of rocket 
exhausts. Measurements based on single particle measurements 
are susceptible to a variety of biasing problems caused by 
correlations between the data rate ana the flow parameters. 
Particle based anemometry also suffers from reduced accuracy 
when used in highly accelerated, high velocity flows. For these 
flows, very small seed particles must be used if the particle 
velocity is to accurately represent the gas velocity. Also, the 
size of the probe volume must be increased so as not to exceed 
the frequency response of the photodetectors. The combination 
of these factors degrades the performance of particle based 
anemometers in a variety of aerospace research programs. This 
has led researchers working with high speed flow to examine 
other approaches to anemometry based on molecular scattering. 


Velocity measurement based on the direct interaction of laser 
light with the gas molecules has several attractive features. 
Both the need for seeding and any particle lag problems are 
eliminated. The simplest of the molecular scattering techniques 
is based on Rayleigh scattering, which is the subject of this 
paper. 

Rayleigh scattering is an elastic process where the internal 
energy of the molecules is not changed. Thus the frequency of 
the scattered light is shifted from the frequency of the incident 
light only by the Doppler effect due to molecular motion. The 
spectrum of the scattered light contains information about the 
velocity distribution function of the gas molecules. In 
particular, the frequency of the peak of the spectrum of the 
scattered light relative to the frequency of the incident light 
represents one component of the gas velocity (in the same 
manner as a conventional reference beam laser anemometer). In 
this paper, we examine the feasibility of using Rayleigh 
scattering for velocity measurements in high speed flows. We 
first look at the fundamental limit on velocity measurement 
accuracy without consideration of the instrument used to 
measure the spectrum of the Rayleigh scattered light. We then 
look at the measurement limits when a Fabry-Perot 
interferometer is used for measuring the spectrum. In this 
analysis we include the effects of spurious laser light, which is 
often a major problem in the application of Rayleigh scattering. 
We analyze both single pass and multiple pass interferometers. 
Multiple pass interferometers offer a much higher frequency 
selectivity and are necessary for situations having a large 
amount of stray laser light, which is not uncommon for many 
aerospace test facilities. 

An example of nitrogen at ambient temperature and pressure is 
presented for a range of Mach numbers and spurious light 
levels. It is shown how the results of this example can be scaled 
for other flow parameters. We then validate our analysis by 
applying it to a Rayleigh scattering diagnostic previously 
developed for use in the exhaust plume of a hydrogen-oxygen 
rocket, where the velocity is about 5000 m/sec. It is shown that 
the results of the analysis are consistent with the experimental 
results. 



Fig. 1 Scattering diagram for Rayleigh scattering example: 
scattering angle 0 S =9O 0 


THEORY 


Consider a linearly polarized plane wave with wave vector k 0 
(fig- J) propagating through a gas with a Maxwellian velocity 
distribution 1 


f ( y ) d3v = i m a3 e~ c V aa d 3 v (1) 

where 


associated with momentum transfer in the scattering, with k s 
and k 0 being the wave vectors of the scattered and incident 
light. The magnitude of K is a function of the scattering angle 
0 S and the wavelength A 


K = I K| = 2k 0 sin( 0 s /2) = ^sin(0 s /2) (5) 


We simplify the equation for the spectrum by introducing 
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where u K is the component of the velocity u in the K direction. 
The spectrum can then be expressed 


S(f) df = 
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We also let 


r . gi p g t,,An [d fi ]. |n , :t ( „ 

be the detected photon arrival rate (i.e., photoelectron count 
rate) of the total Rayleigh scattered light. Here, we used the 
relation IoV sc = P 0 L SC , where Po is the laser beam power and 
L sc is the length of the scattering volume along the beam. The 
overall efficiency factor e accounts for detector quantum 
efficiency and system losses. 


a — (2/tT/m)^ — "most probable speed”, m/sec 
c - v-u 

* = Boltzmann’s constant (1.38x10-23 J/K) 
m = molecular mass, kg 
n = number density, irr 3 
T = gas temperature, K 
u = mean velocity, m/sec 

The power scattered in the k s direction into solid angle dfi in 
frequency interval df is 4 

P s (f)dfdn = I 0 n V s<: [jj^j sin** S(f) dfdO (2) 

and the mean number of photons in time interval At in light 
frequency interval fj to fj+Af is 

< nj > = IonV^Anjdaj sin2;f $(fj) AfAt (3) 

where 

[3^] = Rayleigh differential cross section, m 5 /sr 

I 0 = incident irradiance, W/m 2 

V sc = scattering volume, m 3 

X - an 8 le between incident electric field vector and k s 

S(f) — the normalized spectrum of Rayleigh scattered 

light 

n = solid angle of light collection optics 

For a single component, low density gas (y parameter << 1) 
the spectrum of the scattered light is 2 


Including spurious light at the laser frequency and broadband 
noise, the mean number of counts measured in frequency 
interval f q to f q -f Af in time interval At is 

<n q > =[I\ S(f q )Af 4- r w s qi + r B Af] At (9) 

where I\ is the count rate of detected light at the laser 
frequency (parasitic light scattered from surfaces) and r B is the 
count rate per unit frequency interval of broadband noise (from 
dark noise and background light). The Kronecker delta function 
<5 q i = 1 if q=l (where the subscript 1 refers to the frequency 
interval containing the laser frequency) and = 0 otherwise. 


C ramer- Rao Lower Bound 

In general, the spectrum of the scattered light is a function of 5 
unknown parameters, which we denote by the vector 

a= [r R ,r w ,r B ,u K ,T] (10) 

The measurement can be described by the vector 

n = [n u n 2> ...,n qr ..] (11) 

where the elements of n are the photoelectron counts for equal 
frequency intervals Af in time interval At. The Cramer-Rao 
lower bound for the variance of the estimate of the parameter at 

is 3 


= [r-i]u ( 12 ) 

where T is the Fisher information matrix with elements given by 


S(f)df=^e'( 2 ,rf - K ' u ) 2 /a 2 K 2 df (4) 

The width of this Gaussian function is proportional to the 
square root of the ratio of the gas temperature T to the 
molecular weight m. The entire spectrum is shifted by a 
frequency f D = K-u/2t, where K=k s -k 0 is the wave vector 


fij = - / Jg lnp(n|q) 


da\ da } 


(13) 
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Here, the conditional probability function p(n|a) is the 
likelihood function, which expresses the probability that given a 
set of parameters a, the set of observations n will occur. (The 



maximum likelihood estimate for the parameters a is given bv 
the a that maximizes p(n| a) for a given set of measurements n.) 
For Poisson statistics, which we assume apply to the scattered 
light, the probability of n q counts in one time interval is 

p(n q )= H q e ~ (nq) K) nq ( 14 ) 

where the expected value <n q ) is given by the equation 9. 

The elements of the Fisher information matrix can then be 
written 


q 

If the parameters are uncorrelated, then the variance of the 
estimate of each of the parameters is simply the inverse of the 
corresponding diagonal element of the Fisher matrix, i.e. 


= (srfzjpi = [ 7 ] (2r c R Afj 
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where c s is the speed of sound in the gas and 7 is the ratio of 
specific heats. This shows that, for a fixed Rayleigh scattering 
count, the error in the velocity measurement is proportional to 
the speed of sound in the gas. The relative velocity error is 


Uk 
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where M K = u K /c s can be thought of as the Mach number of the 
velocity component being measured. Thus the relative accuracy 
of the velocity measurement improves as the Mach number 
increases. 


RESULTS 

Example - Nitrogen at ambient conditions 


v„c(&i) = [r-‘]ii = - r L - ( 16 ) 
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For this case, the lower bound for the variance of the velocity 
parameter u' is 



q 


Using equation 9 for <n (1 >, we obtain 



In general, the errors are correlated, and the variance calculated 
from the inversion of the Fisher matrix (eq. 12) should be used. 
The simplified expression given by equation 18 should only be 
used if one is confident that the errors are uncorrelated. 


Lower bound for zero backgroun d and wall scatter 

We can establish simply the lower bound for errors in velocity 
measurement for the case where spurious light at the laser 
frequency and background light are negligible. For this case, 
calculations based on equation 12 show that the velocity error is 
only weakly correlated with the errors in the other parameters. 
Thus the simple form for the Cramer— Rao lower bound for the 
velocity measurement error (eq. 18) can be used, which gives 
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(19) 


We evaluate the sum in equation 19 by allowing Af -* df and 
converting it to an integral. Performing the integration gives 


As an example, consider nitrogen at near— ambient conditions 
(MW = 28 amu, T = 300 K, p = 14.7 psia = 1.014*105 N/m 2 , 
viscosity 4 tj = 1.78*10 -5 N-sec/m 2 , Rayleigh differential 

scattering cross section 5 (dcr/dO) = 7.03* 10 " 32 m 2 /sr). We 
assume a laser wavelength A = 514.5 nm, laser power P 0 = 1 W, 
probe volume length L sc = 1 mm, f/4 collection optics (Q = 0.05 
sr), overall efficiency e = 0.01, \ = 90 * (" s " scattering), and 
scattering angle 0 S = 90* • From the ideal gas law, the molecular 
number density is n = p//cT = 2.45* 10 25 m‘ 3 . Based on these 
values, the detected count rate from equation 8 is P R = 2.23* 10 8 
counts/sec. The speed of sound c 8 = (7*T/m) l/2 = 353 m/sec 
and the most probable molecular speed a = f2^T/m) 1/2 = 422 
m/sec. The scattering wave number K = (49r/A)sin(0 s /2) = 
1.727*107 m-i. 

The y parameter 8 , which indicates the importance of collective 
effects, is 



Although this value of y puts us in the regime where collective 
effects start to become important (and the spectrum is no longer 
strictly Gaussian), we use the Gaussian spectrum for the 
purposes of our example. The results presented here should 
generally be valid for gas densities less than 1 amagat (gas 
density at STP). Note, however, that the deviation from 
Gaussian would be more important if the spectral shape were 
being used to determine the gas temperature. 

Assuming a data acquisition time of At = 1 second, the total 
number of Rayleigh counts is G R = F R At = 2.3* 10 8 , and the 
Cramer-Rao lower bound (for r^r^O) calculated from 
equation 21, is 


<j(uk) = 0.2 m/sec (24) 

Thus, for our example, the theoretically lower bound for the 
velocity error is well within the needs of most experiments. 
Note that this bound does not depend on the velocity, so the 
relative error decreases as the velocity increases. This is in 
marked contrast to particle based anemometry, where the 
relative error increases as the velocity increases. 
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Thus the predicted standard deviation in the measurement of 
the velocity component along K is 


Although the theoretical lower bound for velocity measurement 
is satisfactory, it could only be achieved by using an ideal 
spectrum analyzer. We know that any real instrument used to 
measure the spectrum will give less accurate measurements. In 
the next section we calculate the lower bound for a spectrum 
measured with a particular instrument — the planar mirror 
Fabry— Perot interferometer. 
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Fig. 2 Fabry— Perot interferometer 


M easurement of spectrum using Fa brv-Perot interferometer 

The Fabry-Perot interferometer 7 > 8 (fig. 2) is a simple 
instrument consisting of two parallel, partially transmitting flat 
mirrors. Its transmittance is a function of the frequency of the 
incident light, the angle 0 of the light with the optical axis, the 
refractive index \x of the medium between the mirrors, and the 
mirror separation d. These parameters determine the phase 
change ip between successive reflections, which is (neglecting 
phase change on reflection) 


C0S ° (25) 

For a single ray, the transmittance is given by the Airy function 

lt = FT F sin 2 ( ip/ 2) ( 26 ) 

where F = sin" 2 (2N R /7r), with N R being the reflective finesse. 
The Fabry— Perot can function as a scanning spectrum analyzer 
by measuring the output while varying the mirror spacing d. In 
practice, the light passed by the interferometer is focused with a 
lens called the fringe forming lens. An aperture, with radius r, is 
placed in the focal plane of this lens to define the source area 
from which light is collected. The usable light collection ability 
of a planar mirror Fabry— Perot interferometer is given by the 
etendue U, which is the product of the source area observed and 
the collection solid angle, i.e. 

U = 7rr2 n = TTT 2 *(D/2f L )2 (27) 

where D is the mirror diameter and f L is the focal length of the 
collection lens (for simplicity we assume the focal lengths of the 
collection lens and fringe forming lens are equal). The power 
that enters the interferometer is then given by the product of 
the source radiance and the etendue. 

The spectrum measured with a Fabry— Perot interferometer can 
be written 9 


r.(f d ') = 

s R (f')i*(f d '-f') df' + r w i,(f d ')+ b (28) 

where f'=f-f 0 is the frequency shift relative to the laser 
frequency f 0 . The spectrum is a function of the scanning 
parameter fd /= fo(l~d/d 0 ), which is related to the mirror 
spacing d. The mirror spacing d is the physical scanning 
parameter and varies as a linear function of time from an initial 
spacing d 0 to d 0 +Ad, where Ad is on the order of the 
wavelength of light. The Rayleigh scattered power I\ that 
enters the interferometer has frequency spectrum S R (f'). The 
power F w represents the laser light elastically scattered from 
surfaces. The broadband background light is given by B. 


In equation 28, I R (f<j') and I*(fd') are the instrument response 
functions for Rayleigh scattered light and wall scattered light. 
The instrument function of the Fabry— Perot is the response to a 
monochromatic source with the same spatial distribution as the 
source to be measured 7 . Although I R (f d ') and I w (f d ") are not 
necessarily the same, we assume here, for simplicity, that they 
are equal (i.e., lR(fd / )=fw(fd')=I(fd / ))- We also assume that the 
source has a uniform spatial distribution, so the instrument 
function can be numerically evaluated using 
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i(fd') 

dA 


1 -f F sin 2 


~2 Tf 0 /id { 
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cosfl ip-fd 7 ! 


(29) 


where the integration is over the aperture located in the image 
plane of the fringe forming lens. The angle 6 is the angle of the 
ray in the interferometer measured from the optical axis. As 
shown in figure 2, 9 can be expressed in terms of the position of 
the ray at the aperture (x d ,y d ) 


0 = tan-if- fodj l + ke 

1 l 


(30) 


We form the measured spectrum by setting the mirror spacing 
(and hence the scanning parameter fd') to fixed number of 
spacings that we label q, where q=[ 1,1001. Counts are 
accumulated for a time At at each spacing. The spectrum is 
then in histogram form with the number of counts in the q* h bin 
given by 


<n q > — P m (fq) At 


(31) 


Fabrv-Per ot measurements for nitrogen at NTP 

In this section we analyze the example discussed above for 
nitrogen at near— ambient conditions with the scattering 
geometry shown in figure 1. The predicted error in the 
measured velocity is the Cramer-Rao lower bound given by 
equation 12 (equations 15, 31, 29, and 28 are used in the 
evaluation of equation 12). The integration in equation 29 is 
done using Gaussian quadrature, and the convolution in 
equation 28 is done using fast Fourier transforms. The results of 
this example can be scaled for use with other flows, as is 
discussed in the following section. 

We now consider the measurement of the spectrum using a 
Fabry— Perot interferometer for the above example (nitrogen at 
ambient conditions). The Fabry-Perot is used in a scanning 
mode with the data accumulated in 100 equal time interval bins, 
so At = 0.01 sec with a total data acquisition time of 1 sec! 
Here, we define G R = r R At*(number bins) and G w = 
r w At*(number bins) as the total number of counts over the 
entire data acquisition period. (Note that for the previous 
calculation of the lower bound done without regard to the 
method of measuring the spectrum, we defined G w = r w At, 
where At was the data acquisition period.) We assume a 
Fabry-Perot mirror diameter D = 50 mm and a reflective 
finesse N R = 20. The Fabry-Perot mirror spacing d and 
detector aperture size are free parameters and can be sdected to 
minimize the predicted error in the velocity measurement. For 
a given flow, the uncertainty in the velocity measurement is 
presented in the form of contour plots as a function of free 
spectral range and probe volume diameter. Note that the length 
of the probe volume L sc in equation 8 is determined by the 
projection of the image of the detector aperture onto the beam 
(i.e., L sc is a function of both the scattering angle 9 S and the 
probe volume diameter). We assume that the probe volume 
diameter is equal to both the diameter of the image of the 
circular detector aperture and the diameter of the laser beam. 
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NITROGEN 
T = 300 K 
P = 14.7 psia 
Gw/Gr = 1 
Mach 1 


LEVEL ERROR 
8 6.00 m/sec 

7 5.50 

6 5.00 

5 4.50 

4 4,00 

3 3.50 

2 3.00 

1 2.50 
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Fig. 3 Predicted velocity measurement error for single pass 
Fabry-Perot interferometer as a function of probe 
volume diameter and free spectral range 



Each contour plot is based on a 10*10 grid of computed values. 
The computation time for each contour plot, using a 33 MHz 
386 personal computer, was about 700 seconds. 

Results were obtained for four Mach numbers (0.2, 1, 2, 4 ), six 
values of the ratio of wall scatter to Rayleigh scatter ratio 
G w /G r (0.1, 1, 10, 100, 1000, and 10000), and 1, 3, and 5 pass 
Fabry-Perot configurations. (The multi-pass transmittance is 
given by the single pass transmittance, eq. 26, raised to the N th 
power, where N is the number of passes 10 .) Figure 3 shows a 
representative contour plot for Mach 1, G W /G R = 1, and single 
pass. The minimum values of the velocity uncertainty from 
these contour plots are plotted as a function of G*/G R in figure 
4. Note that the minima are rather broad in figure 3, so the 
exact values of mirror spacing and detector size are not critical. 
(However, the error does increase more rapidly for free spectral 
range smaller than optimum.) Even for a small amount of wall 
scatter, the predicted uncertainty in the velocity measurement is 
about 2 m/sec. This is approximately an order of magnitude 
larger than the lower bound calculated for measurements 
obtained with an ideal spectrum analyzer. The primary cause of 
this difference is that the scanning Fabry-Perot interferometer 
uses only a fraction of the incident light at any time (that light 
whose frequency falls with the passband of the interferometer), 
whereas the ideal spectrum was assumed to use all of the 
incident light for the entire data acquisition period. As one 
would expect, as the amount of wall scattered light increases, 
the accuracy decreases. 

Figure 4a is for a single pass Fabry-Perot and shows a large 
degradation in accuracy for G W /G R > 10. Figures 4b and c show 
similar data for 3 pass and 5 pass configurations. These show 
the great improvement in accuracy achievable using multi-pass 
configurations in experiments with large amounts of wall 
scattered light. It should be noted that we assumed the same 
Fabry-Perot mirror diameter for these multi-pass 
configurations as was used for single pass. Although this could 
be achieved, in principle, by operating interferometers in series, 
the usual technique is to use a single interferometer with 
retroreflectors to obtain multiple passes 10 . This method greatly 
reduces the effective mirror aperture, which decreases the 
amount of usable Rayleigh scattered light. In addition, the 
losses in the Fabry-Perot increase with the number of passes. 
For example, a single pass interferometer with 75% peak 
transmission would only have 42% and 24% transmission when 
operated in 3 pass and 5 pass configurations, respectively. 
Finally, for sake of comparison, the predicted velocity 




Fig. 4 Predicted minimum velocity measurement error as 
function of ratio of wall scatter light to Rayleigh 
scatter light for nitrogen at 300 K: (a) single pass 
Fabry-Perot (b) 3 pass (c) 5 pass 


5 



NITROGEN 
T = 300 K 
P = 14.7 psia 



> 


0 00 L — L. 

0 25 


0.50 0.75 1.00 

PROBE VOLUME DIAMETER, mm 


Fig. 5 Predicted uncertainty in velocity measurement for 
ideal spectrum analyzer 


measurement error for an ideal spectrum analyzer is shown in 
figure 5 (with the error plotted as a function of probe volume 
diameter). Note that an ideal spectrum analyzer would not be 
affected by wall scatter light that is located in an infinitely 
narrow bandwidth. 


In general, the predicted error in the velocity measurement is 
proportional to the speed of sound in the gas and inversely 
proportional to the square root of the number of detected 
photons. This is the relation found for the case of the 
theoretical lower bound (eq. 21), and was verified for 
calculations made for the Fabry— Perot with large amounts of 
wall scatter and for multi-pass configurations. This relation 
permits the results obtained for the above example to be 
extended to other situations. We consider two examples. 

The first example of scaling is a low density (0.01 amagat), 
Mach 4, nitrogen flow having a static temperature of 900 K. 
Compared to the ambient condition example, the Rayleigh 
counts would be decreased by a factor of 100, and the ratio 
W/G r would increase by the same factor (assuming the same 
amount wall scatter), say from 1 to 100, and the sonic velocity 
would increase by a factor of 3 1/3 due to the higher temperature. 
tor a single pass Fabry-Perot, using figure 4a, we see that the 
velocity error would be (3.3)(10)3" 3 = 57 m/sec, or about 2.4%. 

The second example of scaling is a Mach 1, nitrogen flow at 
ambient temperature and pressure in an experiment with a large 
amount of wall scatter, such as one would typically have in an 
internal flow. Because of the assumption of a large amount of 
wall scatter, we select a 5 pass Fabry-Perot. We assume that 
the effective mirror diameter is 7 mm, and the losses in the 
rabry— Perot are such that the efficiency is only 0.25% instead 
of the 1% in the above example. These differences cause the 
Rayleigh count rate to be reduced by a factor of about 200 
which corresponds to an increase of the predicted error in the 
velocity by a factor of 14. Using this factor with figure 4c, the 
predicted error is about (3)(14) = 42 m/sec, or 12%. Here, we 
might consider extending the data acquisition time (e.g., 
increasing the time from 1 second to 10 seconds would decrease 
the predicted error to 13 m/sec, or 4%). 



Fig. 6 Scattering diagram for hydrogen— oxygen rocket 
measurements; scattering angle 0 S =152.8* 
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T = 1 000 K 
P * 0.08 psia 
Gw/Gr = 1 
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Fig. 7, Predicted velocity measurement error for 
hydrogen— oxygen rocket exhaust measurements 



Comparison with experimental results from Ht-O, rocket 

We now apply our prediction analysis to a Rayleigh scattering 
diagnostic developed for a small (100 N thrust) 
hydrogen-oxygen rocket test facility. This system and the 
experimental results obtained with it are described in reference 
11. The relevant parameters for the H 2 0 exhaust gas are: 
MW=18, 7=1.25. p=0.08 psia, T=1000 K, V=5000 m/sec, 
05=152.8-, (d(r/dn)=5.01xl0‘ 33 m 3 /sr. The single pass 
Fabry-Perot used had an entrance aperture of 64 mm and a 
mirror spacing of 7.0 mm (FSR = 21.4 GHz). A 1 W argon-ion 
laser was used at A = 514.5 nm. The light was linearly 
polarized perpendicular to the scattering plane, so *=90’ . As 
shown in the scattering diagram (fig. 6), the measured velocity 
component was about 75‘ from the actual velocity, so the Mach 
number of the measured component of the flow was M* = 1.7. 
The probe volume diameter was 2.2 mm and length was 4 mm. 
Spectral data were collected for 90 seconds. Assuming an 
overall efficiency of 1%, the expected detected Rayleigh 
scattering count rate (from eq. 8) is T R = 3300 counts/sec. For 
measurements taken on the centerline of the exhaust plume, the 
measured Rayleigh count rate was about 5000 counts/sec. The 
background count rate was about 250 counts/sec, ana the ratio 
Gw/G r typically was about 1. 

Applying our analysis to this flow and Fabry-Perot 
configuration gives the predicted error shown in figure 7. From 
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this figure, the optimum probe volume diameter is 2.6 mm 
(close to the value of 2.2 mm used in the experiment). From 
figure 7, we see that the optimum free spectral range is 14 GHz, 
which is somewhat less than the experiment value of 21.4 GHz 
However, this would only reduce the error by about 20%. The 
predicted Rayleigh count rate is about 4800 counts/sec, which is 
close to the value obtained in the actual measurements. This 
indicates that our assumption of 1% overall efficiency is 
reasonable. 


CONCLUDING REMARKS 

The prediction analysis applied to the measurement of velocity 
based on the spectral measurement of Rayleigh scattered light 
with a planar mirror Fabry-Perot interferometer indicates its 
feasibility for a variety of flow conditions. In experiments where 
the amount of wall scattered light cannot be reduced to low 
levels, the use of a multi-pass Fabry-Perot offers a means to 
achieve adequate frequency selectivity. 

The use of a scanning Fabry-Perot interferometer precludes 
instantaneous velocity measurements. Even with relatively high 
Rayleigh scattering levels, several seconds per measurement 
point will probably be required. In rotating systems, the 
measurements would need to be correlated with the angular 
position of the rotating part, thus further increasing the data 
acquisition time. 

The material presented here considered only the measurement of 
a single velocity component. The method, however, can be 
extended to multiple components in a straightforward manner 
by adding either additional receiving optics, or by adding 
additional input beams. 

It is important to note that the velocity measurements described 
in this paper are not sensitive to the exact molecular 
composition of the gas (because the mean velocity of the gas is 
determined by the frequency of the spectral peak of the Rayleigh 
scattered light). Furthermore, if a large, uniform concentration 
of small particles is present in the flow, the velocity 
measurement can be made using the light scattered from the 
particles. (Examples of particle scattered light used to measure 
velocity using a confocal Fabry-Perot are presented in 
references 12 and 13). If a small number of relatively large 
particles are in the flow, they will scatter large amounts of light 
compared to the molecular Rayleigh scattering, and the spectral 
data can be screened to eliminate the anomalously large values 
caused by the particles. 

Finally, we would like to emphasize that: (1) since the predicted 
relative error decreases with Mach number, the Rayleigh 
scattering technique is particularly suited to supersonic and 
hypersonic flow studies; and (2) because the velocity measured is 
the actual molecular velocity, highly accelerated flows, including 
flows with shocks, can be accurately mapped. 
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